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The nature of cosmological dark matter remains mysterious. Recently, Spergel and Steinhardt [1] 
- have revived suggestions [2, 3] that dark matter may be strongly self-interacting, i.e., collisional, to make 

■ galactic halos less dense. We show here an important by-product: for reasonable values of particle mass 
; j ' and collisional cross-section, galaxy cores would quite naturally grow, within them, massive black holes 

■ 1O^-1O^M0, having the scaling observed by Magorrian et al. [4], yisH « ^^'^ai- 



The physical picture is quite simple. Assuming a normal power spectrum of perturbations, dark matter 



04 ' halos begin to form in earnest in the redshift range Z = 30 — s- 20, with star formation commencing in a 

> 

oo 
in 

psj ■ massive stars formed at these early epochs (Abel et al. [7]) will have several dramatic effects on subsequent 



significant way during the interval Z = 20 ^ 10 (Ostriker and Gncdin [5], Haiman, Recs and Loeb [6]). The 



cosmic evolution: they emit UV radiation copiously, which begins to reheat and reionize the intergalactic 

' medium; they ultimately explode, contaminating their environment with the first heavy elements; and, 

Q ' most importantly for our present purposes, their cores implode to leave black holes (Arnett [8]), having 

■ masses ^1/4 of the original stellar mass. The existence of a population 111 of high mass stars which 

Cd ■ 

would have left black hole remnants is attested to by the non-zero {Z/Zq k, 10^"^) floor to the metallicity 
><'... 

; I ■ distribution observed ubiquitously (9) indicating very early and widespread contamination of the universe 

, 

by ejecta from high mass supernovae. How such stellar mass black holes grow to the supermassive size 
(10^ < Mbh/Mq < 10^) seen in galactic nuclei [4] remains unknown, although rather general considerations 
(Rees [10]) indicate the plausibility of their formation. Furthermore, the quantitative scaling observed [4] 
Mbh oc Mgai oc Vj,'* ^^'*'^ (where Vc is the galactic rotation velocity) is an outstanding mystery. 

The 25 Mq black hole remnant from a 100 Mq {H, H^) star will immediately begin to accrete collisional 
dark matter in the core of the dark halo within which it forms. To frame the discussion, consider the 
simplest case first. At high densities, the dark matter will behave like an adiabatic gas and accrete as per 
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the classic work of Bondi [11]: 



Mbh = Mac = 47rriC^pA with ta = GMbh/C\, 



(1) 



where {pA, Ca) are the density and sound speed in the ambient dark matter fluid. 

Accretion is treated in the subsequent discussion as if it were quasi-spherical, quite different from the 
normal treatment of accretion from a rotationally flattened disc. The reason for this is that for normal 
baryonic matter, electromagnetic radiative losses are usually eflicient enough so that energy loss dominates 
over angular momentum loss in environments of accretion, which leaves the contracting matter in a disc 
whose further evolution is limited by the rate at which viscosity, feeding from the quasi-Keplerian differential 
rotation, can transport energy and angular momentum outwards [12]. But in the non-radiative case, both 
energy and angular momentum can be transported only by collisions (gravitational or physical). Then, as 
shown by Goodman [13] and others, both transport processes occur on approximately the coUisional time 
scale {i.e., the viscous and conductive time scales are comparable), with angular momentum transported 
outwards rapidly enough so that the central regions remain quasi-spherical. 

Stone et al. [14] have investigated a similar situation of non-radiative accretion in the limit of significant 
rotation, low viscosity and negligible conduction. They convincingly find both convection and outflow to be 
important, but it is not likely that their solution applies to the high viscosity, high conductivity, low rotation 
case considered here. In any case, most of the accretion occurs in the second, optically thin regime which 
we will address shortly. 

We thus return to the simple illustrative example of Bondi - accretion in a quasi-spherical dark matter 
halo. The core region of which this matter is composed is likely isothermal in its prior structure due to both 
collisional and violent (Lynden-Bell [15]) relaxation and thus has a profile: 



(We consider a more general profile subsequently.) 

Integrating equations (1) and (2), one finds that the central massive black hole grows at the speed of 



PA(r) = C\l (27rGr2) . 



(2) 
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sound to reach, in time t, a mass GMbh (t) = 'iC\t, where we note in passing that it will vacuum up all the 
local baryonic components as well as the core dark matter. 

This phase of rapid growth can only persist until the accretion radius has grown to the point far enough 
from the center whereby the dark matter mean free path approaches the accretion radius. Then a transition 
will occur to slower, diffusively limited growth. At this time, ti, mean free path A = mp/{pt^(T) = ta = 
GMbhmIC% or 



where {ap,mp) are the dark matter self-interaction scattering cross-section and particle mass respectively. 
At this time, the core radius is r^^ = crj,C^/(27rGmj,). An alternate way of writing equation (3) in terms of 
dimensionless numbers is instructive: 



rUe 2 \Gml) \mp) \crT ) \ c 



(4) 



where the first quantity in parenthesis is the classic "big number" representing the ratio of electromagnetic 
to gravitational forces (~ 10''^), or is the electromagnetic Thompson cross-section and c is the speed of 
light. To translate this into physical units, we take rUp = 100GeV,a/aT = 10~^ and Ca = WOkm/s, then 
Mbh = 3 X IO^Mq, an interesting value. 

The next phase of slower growth has been treated by several authors. A cusp forms about the black hole 
approximately described (for gravitational interactions) by the classic Bahcall-Wolf [16] solution, as modified 
by loss cone effects (Ostriker and Tremaine [17]) in the inner parts. The rate of accretion is determined by 
the rate at which particles are scattered into the loss cone, which begins at the Bondi accretion radius. Thus, 
it is initially the Bondi accretion rate multiplied by the probability of strong scattering for a particle orbiting 
at the Bondi radius: P = (jC\/ {GnipGMBH), to give an accretion rate GMbh = 2aC\/{GMBH){Gmp). 
This produces a mass growing, after the transition, as the square root of the time, giving for the present 
time tH- GMbh = \^^o-C\tH / {GrUp) , which gives Mbh = 4.1 x 10^ for the previously quoted parameters, 
and tn = 10^° yrs. This produces a somewhat too large value for Mbh, but one with approximately the 
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correct scaling on Ca- In fact, this solution cannot be extrapolated to late times or to very small values 
of a/mp, because at some point accretion onto the black hole will be limited by the mean free path in the 
vicinity of the Schwarzschild radius, after which accretion only occurs within a loss cone. Then the estimate 
is reduced by a factor of (C^/c) to give (Ostriker and Tremaine [17]) 



Since the fundamental particle physics is quite uncertain in any case, we can best parameterize our 
ignorance by defining r] = a/Gnip and, noting that the requirement of Spergel and Steinhardt [1] that the halo 
be optically thick to collisions at a radius ri is equivalent to 1 = Annicj = 2C\(j/'KGmpr\, or 77 = 0.57rri/C^ 
(it is trivial to generalize these definitions to the likely case that the cross-section is velocity dependent), 
we find GMbh = 2.Qy/C\tHri/c'^ or Mbh/Mq = 4.9 x 10^c[^o^^,io^ifepc- The range for n quoted by 
Spergel and Steinhardt is approximately O.lkpcCiQQ <ri < O.lMpcCiQQ, or 0A5cmP/g < a /nip < ibOcm^/g. 
Then, substitution into equation (5) gives Mbh = 7.06X 106(CT/mp)i/2v;^/2^i^/2g solar masses or, taking the 
midpoint of the proposed range a/rUp = lA.2cm? /g^ the result Mbh = 2.7 x 10'^yg^(oo*H,i5-^©' with a range 
about this (10^'^ ''^) encouraged by the Spergel Steinhardt analysis. We replace the speed of sound Ca with 
the circular velocity V'^ = 2C\ in the above relation to connect more conveniently to normal astronomical 
measurements. For both our own galaxy where K,ioo ~ 2, and for M87 for which Vc,wo « 7, the resulting 
black hole masses are too large, as compared to observed estimates [4] by a factor of several hundred. This 
should be treated as a remarkable agreement given the crudeness both of the analysis and of the estimated 
collisional DM properties. Accretion in the optically thin limit is only significant for the lower part of this 
mass range. 

Clearly, the assumption of a single stellar mass black hole in the central region of the galaxy is an over 
simplification, since multiple seeds are expected. Given the instability of the initial growth, one BH is likely 
to dominate and will eat or eject the others [18]. This complication is not likely to significantly alter the 
estimates of the final black hole mass. 




(5) 



where the numerical constant / is approximately / = [5 + 21n(c/CA)]/27r w 2.60. 
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It is reasonable to ask how the solution changes if we relax the assumption that the profile is that of a 
singular isothermal sphere for which the density profile (eq. 2) is p oc r~" with a = 2. The widely adopted 
NFW [19] profile for dark matter halos takes a = 1 in the inner parts, and other authors find typical values 
near a = 3/2 (see [20, 21] for references). One can show [17] that, in the more general case, if the galactic 



M oc 1/''+'', then the final black hole mass will scale with the galaxy mass as Mbh AfJ^; ) where 



Remarkably enough, for the case of greatest interest, where p + (5 = 9/2, then K = 1 regardless of the value of 
OL. Thus, the conclusion reached in the case of a cusp which is initially that of a singular isothermal sphere 
(a = 2), that Mbh oc Mgai-, the observational result found by Magorrian et al. [4], is hkely to be very close 
to that generally obtained for values of the initial cusp parameter departing moderately from this value. 

In the most plausible [20, 21] range for the parameter a, [a. = 1.3±0.2), the core profile would evolve due 
to the collisional effects, initially reducing the central density [1, 22] and lowering the rate of growth of the 
central BH, i.e., the evolution of the core must be treated simultaneously with the growth of the central BH. 
Quinlan's [22] numerical work indicates that the collisional evolution of the core ultimately approaches the 
singular core collapse {a — 2.25) solution [23] which would, if it were reached, greatly accelerate the growth 
of the central black hole. However, for the acceptable range of the parameter a/m, this state would not 
be reached within a Hubble time for all except the most extreme systems (see also Burkert [24]). Also, the 
rotation of the dark matter halo will not be trivial, but, since the viscous and angular momentum (outward) 
transport times are the same as the other relevant time scales [13], this is not likely to provide a significant 
barrier to accretion. Both of these effects reduce the estimated mass of the BH to a level below that given 
by equation [5]. Nevertheless, the limit on the observed sizes of central galactic black holes (perhaps below 
the Magorian et. al [4] estimates) probably restricts a /nip to the lower end of the range proposed by Spergel 
and Steinhardt: a /nip < 1cm? /gm. However, only detailed numerical calculations (now in progress) will 
be able to establish a more precise bound. 



luminosity scales as L a V^, the mass-to-light ratio scales as (M/L) cx , 



so that the mass scales as 




(6) 
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It is quite possible that physical processes not included in this elementary treatment could substantially 
inhibit the growth of black holes in the Spergel-Steinhardt scenario, the required values being somewhat 
smaller than given by equation [5]. We must conclude that the existence of strongly self- interacting dark 
matter has the exciting potential for leading to the growth of central massive black holes in normal galaxies 
with the observed scaling parameters. 

Three corollary consequences should be noted. First, since the hypothesized dark matter particles do 
not have radiative interactions, accretion of them will not produce an electromagnetic luminosity output, 
Lem, breaking the assumed link J Lemdt = e,;,n^M/c'^, with normal estimates of eem ~ 0.1. Alternatively 
phrased, if most of the accreted matter is dark matter, then a low efficiency (eem << 1) is to be expected. 
Second, since the mean free path of the particles is comparable to the system size, the dark matter fluid 
will be extremely viscous, with dynamical consequences that may be imagined. Finally, dark matter 
galactic halos in clusters of galaxies will tend to evaporate due to heat transfer from the hotter, cluster 
dark matter. Preliminary estimates of the significance of this effect again limit one to consider the range 
(a/rUp) < Icrri^ / gm. 
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